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Using  a  specially  designed  current  distribution  measurement  gasket  in  anode  and  thin  thermocouples 
between  the  catalyst  layer  and  gas  diffusion  layer  (GDL)  in  cathode,  in-plane  current  and  temperature 
distributions  in  a  proton  exchange  membrane  fuel  cell  (PEMFC)  have  been  simultaneously  measured. 
Such  simultaneous  measurements  are  realized  in  a  commercially  available  experimental  PEMFC.  Exper¬ 
iments  have  been  conducted  under  different  air  flow  rates,  different  hydrogen  flow  rates  and  different 
operating  voltages,  and  measurement  results  show  that  there  is  a  very  good  correlation  between  local 
temperature  rise  and  local  current  density.  Such  correlations  can  be  explained  and  agree  well  with  basic 
thermodynamic  analysis.  Measurement  results  also  show  that  significant  difference  exists  between  the 
temperatures  at  cathode  catalyst  layer/GDL  interface  and  that  in  the  center  of  cathode  endplate,  which 
is  often  taken  as  the  cell  operating  temperature.  Compared  with  separate  measurement  of  local  cur¬ 
rent  density  or  temperature,  simultaneous  measurements  of  both  can  reveal  additional  information  on 
reaction  irreversibility  and  various  transport  phenomena  in  fuel  cells. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

It  is  well  known  that  current  density  is  not  uniform  in  a  proton 
exchange  membrane  fuel  cell  (PEMFC),  and  such  non-uniformity 
has  significant  negative  effects  on  both  its  overall  performance 
and  durability.  Thus,  in  situ  measurements  of  local  parameters 
are  essential  for  better  understanding  of  this  issue  and  for  opti¬ 
mization  of  PEMFC  design  and  operation.  Various  in  situ  local 
diagnostic  techniques  have  been  developed,  including  measure¬ 
ments  of  current  distribution,  temperature  distribution,  membrane 
water  content  distribution,  liquid  water  distribution  and  species 
concentration  distribution,  etc.  [1]. 

Among  the  various  in  situ  diagnostic  techniques,  much  attention 
has  been  paid  to  local  current  measurements  since  local  current 
density  is  a  straightforward  indicator  of  local  performance.  Other 
valuable  information,  such  as  reactant  shortages,  membrane  dry¬ 
ness  and  electrode  flooding,  can  also  be  deducted  approximately 
from  the  distribution  of  local  current  density.  Segmentation  [2-14] 
or  semi-segmentation  [15-17]  of  current  collector/flow  field  plate 
and  fabrication  of  subcells  electrically  isolated  from  main  cell 
[18-21]  are  widely  used  approaches  to  current  distribution  mea¬ 
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surement  in  PEMFC.  Adoption  of  embedded  resistors  [18,22]  or 
Hall  sensors  [23]  were  also  proved  effective.  Indirect  deduction 
from  current-induced  magnetic  field  measurement  [24,25]  or  from 
potential  distribution  measurement  [26,27]  and  separate  measure¬ 
ment  of  current  density  under  the  channel  and  the  shoulder  with 
partially  catalyzed  electrode  [28,29]  are  recently  developed  meth¬ 
ods  with  special  features. 

A  number  of  temperature  distribution  measurement  methods 
have  been  also  developed.  Insertion  of  commercially  available 
or  specially  fabricated  fine  temperature  sensors  into  different 
positions  in  PEMFC,  including  thermocouples  [30-38],  thermistor 
[39,40],  in-fiber  Bragg  grating  sensor  [41]  and  bandgap  tem¬ 
perature  sensor  [42],  is  a  most  widely  used  approach.  Infrared 
thermography  is  another  important  technique  for  temperature  dis¬ 
tribution  measurement  [43-46].  Measurement  of  gas  temperature 
along  the  flow  channel  with  tunable  diode  laser  absorption  spec¬ 
troscopy  was  also  demonstrated  [47]. 

It  is  well  known  that  local  temperature  in  PEMFC  is  intimately 
related  to  local  performance  because  temperature  has  significant 
effects  on  electrochemical  reaction  kinetics  while  local  current  den¬ 
sity  determines  the  local  heat  generation  rate  and  thus  the  local 
temperature.  Furthermore,  some  critical  phenomena  in  PEMFC, 
such  as  water  flooding,  membrane  dehydration,  and  even  cell  fail¬ 
ure  due  to  overheating,  all  depend  on  both  local  temperature  and 
local  current  density.  Therefore,  compared  with  independent  mea¬ 
surement  of  either  local  current  distribution  or  local  temperature 
distribution,  simultaneous  measurement  of  both  can  provide  addi- 
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tional  information  and  lead  to  better  understanding  the  various 
phenomena  in  a  fuel  cell.  Besides,  the  relationship  between  local 
temperature  and  current  density  obtained  from  the  simultane¬ 
ous  measurement  can  reveal  more  fundamental  information  about 
reaction  irreversibility  and  various  transport  phenomena  in  a  fuel 
cell. 

Although  numerous  in  situ  measurements  of  current  distri¬ 
bution  and  temperature  distribution  [2-47]  have  been  reported, 
simultaneous  measurements  of  current  and  temperature  distribu¬ 
tion  in  PEMFC  are  very  limited.  Hakenjos  et  al.  [43,44]  designed 
a  special  PEMFC  allowing  simultaneous  evaluation  of  current, 
temperature  and  water  distributions  in  the  fuel  cell.  Current  dis¬ 
tribution  was  measured  with  segmented  anode  flow  field  and 
temperature  was  measured  with  infrared  thermography  through 
an  optical  window  in  the  cathode  end  plate.  Due  to  the  limita¬ 
tion  of  the  technique,  it  was  only  possible  to  measure  temperature 
distribution  on  the  surface  of  cathode  gas  diffusion  layer  (GDL). 
Maranzana  et  al.  [35]  developed  a  single  PEMFC  to  measure  cur¬ 
rent  and  temperature  distribution  along  the  gas  channels.  They 
observed  that  the  temperature  field  was  strongly  correlated  with 
the  current  distribution  and  both  depended  strongly  on  water  man¬ 
agement.  However,  the  thermocouples  were  inserted  in  the  ribs  of 
the  flow  field  plates  and  thus  it  was  not  possible  to  measure  the 
temperature  close  to  the  catalyst  layer. 

A  large  portion  of  the  heat  generated  in  a  PEMFC  occurs  in 
the  cathode  catalyst  layer  and  its  local  temperature  is  inherently 
related  to  its  local  current  density  and  overpotential.  Thus  in  order 
to  understand  the  relationship  between  local  temperature  and  cur¬ 
rent  density,  for  simultaneous  measurements  of  temperature  and 
current  density  the  ideal  location  for  thermocouples  is  the  cath¬ 
ode  catalyst  layer.  Therefore,  the  objectives  of  this  work  are  (a)  to 
develop  a  technique  to  simultaneously  measure  local  current  den¬ 
sity  and  the  temperature  in  the  cathode  catalyst  layer  and  (b)  to  find 
the  relationship  between  local  temperature  and  current  density. 

2.  Experimental  methodology 


Table  1 

Geometric  parameters  of  the  experimental  PEMFC. 


Channel  length  (mm) 

39 

Channel  width  (mm) 

2 

Channel  depth  (mm) 

1 

Shoulder  width  (mm) 

2 

Flow  field  plate  thickness  (mm) 

13 

Current  collector  thickness  (mm) 

1.5 

Endplate  thickness  (mm) 

19 

Number  of  channels 

10 

Number  of  shoulders 

9 

MEA  active  area  (cm2) 

4x4 

Membrane  thickness  (mm) 

0.05 

CL  thickness  (mm) 

0.01 

GDL  thickness  (mm) 

0.20 

China)  and  two  identical  carbon  paper  GDLs  with  microporous  layer 
(Spectracarb™  2050-L  carbon  paper,  Spectracorp,  USA).  Thick¬ 
ness  of  the  membrane  is  about  50  p,m  and  that  of  GDL  is  200  p,m. 
Catalyst  loading  of  the  CCM  is  0.3  mg  Pt  cm-2  at  the  anode  and 
0.5  mgPtcnrr2  at  the  cathode.  Geometric  parameters  of  the  MEAs 
are  shown  in  Table  1. 

A  pair  of  graphite  flow  field  plates  with  serpentine  flow  field 
pattern  was  used  in  the  fuel  cell.  The  flow  field  at  the  cathode  is 
identical  to  that  at  the  anode  and  parallel  flow  (co-flow)  is  used. 
Channel  width,  channel  depth  and  shoulder  width  are  2  mm,  1  mm 
and  2  mm,  respectively.  There  are  10  channels  and  9  shoulders  in 
every  flow  field  plate.  Details  of  the  geometric  parameters  of  the 
flow  field  plates  are  shown  in  Table  1. 

A  fuel  cell  test  station  (FCTS-1 6,  Fuel  Cell  Technologies,  Inc.,  USA) 
is  used  to  control  operating  parameters  of  the  fuel  cell,  including 
gas  flow  rates,  gas  humidification  temperatures,  backpressure,  cell 
operating  temperature  and  cell  operating  voltage.  Details  about  this 
fuel  cell  test  station  are  described  previously  [17]. 

Hydrogen  (purity  >99.99%)  is  supplied  as  fuel  by  a  hydrogen  gen¬ 
erator  from  pure  water  electrolysis  and  air  is  supplied  as  oxidants 
by  an  air  pump  with  filtering  and  drying  system. 


2.1.  The  experimental  system 


2.2.  Measurement  method 


A  single  PEMFC  with  an  active  area  of  4  cm  x  4  cm  was  used 
for  all  measurements  in  this  work.  Fig.  1  shows  schematically 
the  PEMFC  assembly.  MEA  in  the  fuel  cell  consists  of  a  catalyst 
coated  membrane  (CCM,  Shanghai  Pearl  Hydrogen  Power  Source, 
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Temperature  distributions  in  the  experimental  PEMFC  are  mea¬ 
sured  by  10  very  thin  T  type  thermocouples  (RQB-0.08,  Captec 
Enterprise,  France)  placed  between  catalyst  layer  and  GDL  in  cath¬ 
ode  side  considering  that  most  heat  is  generated  in  the  cathode 
catalyst  layer  which  should  have  the  highest  temperature  in  the 
through-plane  direction.  Together  with  Fig.  1,  Fig.  2  shows  the 
positions  of  the  10  thermocouples.  It  can  be  seen  that  there  is 


Fig.  1.  Schematic  diagram  of  the  experimental  PEMFC. 


Fig.  2.  Schematic  positions  of  thermocouples. 
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Fig.  3.  Photograph  of  the  current  distribution  measurement  gasket  placed  on  the 
anode  flow  field  plate. 


a  thermocouple  above  each  channel  and  the  measuring  tip  is  at 
the  center  of  the  channel  along  flow  direction.  The  thermocouples, 
insulated  with  perfluoroalkoxy  (PFA),  having  a  cross-dimension  of 
0.1  mm  (thickness)  and  0.2  mm  (width)  including  insulation  layer, 
are  very  thin  compared  with  the  channel  width  of  2  mm.  Temper¬ 
atures  measured  by  the  10  thermocouples  are  recorded  by  a  16-bit 
multi-channel  data  acquisition  device  (MX100,  Yokogawa  Electric 
Corp.,  Japan)  at  a  sampling  rate  of  one  time  per  second  for  each  local 
temperature. 

In  many  studies,  the  temperature  measured  by  a  thermocouple 
placed  in  the  cathode  endplate  (CEP),  as  shown  in  Fig.  1,  is  taken 
as  the  cell  operating  temperature.  In  this  study,  this  temperature  is 
named  cathode  endplate  temperature  (  TCep )  to  differentiate  it  from 
the  local  temperatures  measured  at  the  cathode  catalyst  layer/GDL 
interface. 

To  provide  specific  boundary  conditions  for  the  temperature 
measurement,  cathode  endplate  temperature  of  the  cell  is  kept 
stable  with  a  temperature  controller  (Series  1600,  Love  Controls 
Division,  Dwyer  Instruments,  Inc.,  USA),  which  controls  two  heating 
rods  inserted  in  the  endplates  for  heating  and  a  fan  for  cooling. 

Current  distributions  in  the  experimental  PEMFC  are  mea¬ 
sured  with  the  current  distribution  measurement  gasket  technique, 
details  of  which  are  described  in  previous  publication  [17].  Fig.  1 
shows  briefly  how  the  technique  works  and  Fig.  3  is  a  photograph  of 
the  current  distribution  measurement  gasket  placed  on  the  anode 
flow  field  plate.  As  is  shown  in  Figs.  1  and  3,  the  current  distribution 
measurement  gasket  fabricated  with  printed  circuit  board  (PCB) 
technology  is  placed  between  anode  GDL  and  anode  flow  field  plate, 
with  current  conducting  gold  plated  copper  strips  in  direct  con¬ 
tact  with  the  GDL.  The  current  conducting  strips  are  insulated  from 
each  other  and  from  the  flow  field  plate  because  the  gasket  sub¬ 
strate  is  nonconductive.  Therefore,  instead  of  passing  through  flow 
field  plate  in  a  whole,  the  individual  current  in  the  anode  passes 
through  each  conducting  strip  separately.  Thicknesses  of  the  gas¬ 
ket  substrate  layer,  the  copper  layer  and  the  gold  layer  are  about 
0.2  mm,  0.03  mm  and  0.002  mm,  respectively. 

Different  from  the  experimental  system  in  previous  work  on 
current  distribution  measurement  [17,48-50],  in  which  a  multi¬ 
channel  potentiostat  was  used  to  control  the  operating  voltages 
and  measure  the  local  currents,  in  this  study  the  electronic  load 
(N3306A,  Agilent,  USA)  embedded  in  the  fuel  cell  test  station  is 


used  to  control  the  operating  voltage  and  Hall  sensors  (LTSR  6- 
NP,  LEM,  Switzerland)  are  used  to  measure  the  local  currents  in  a 
non-contact  way.  The  range  of  the  Hall  sensors  is  set  to  be  2  A  (low¬ 
est  available  range)  for  high  accuracy.  The  voltage  output  of  the 
Hall  sensor  is  proportional  to  the  current  passing  through  it  and 
it  is  converted  to  current  value  by  the  data  acquisition  software. 
Both  voltage  and  converted  current  values  of  each  Hall  sensor  are 
recorded  by  the  same  data  acquisition  device  used  for  temperature 
recording  at  a  sampling  rate  of  twice  per  second,  thus  simulta¬ 
neous  measurements  of  current  density  and  temperature  were 
ensured. 

2.3.  Experimental  data  reduction 

As  can  be  seen  in  Figs.  2  and  3,  there  are  10  gas  channels  and 
9  shoulders  in  the  flow  field.  For  all  experiments  in  this  work,  the 
experimental  PEMFC  stands  vertically  with  reactant  gas  inlets  on 
the  top  and  outlets  on  the  bottom  of  the  flow  field.  As  is  shown  in 
Fig.  2,  the  10  thermocouples  are  numbered  from  T/C  1  to  T/C  10 
along  the  direction  from  inlet  to  outlet,  corresponding  to  the  10 
channels.  Fig.  3  shows  that  the  1 1  current  conducting  strips  in  the 
current  distribution  measurement  gasket  are  numbered  from  1  to 
1 1 ,  with  number  2  to  number  1 0  corresponding  to  the  9  shoulders. 

Among  the  1 1  current  conducting  strips  shown  in  Fig.  3,  each  of 
the  9  strips  from  number  2  to  number  10  collects  current  from  one 
shoulder  and  two  half-channels,  which  represents  an  active  area 
of  1.6  cm2,  and  the  local  current  density  can  be  calculated  from 
measured  local  current  divided  by  the  area  of  1.6  cm2.  For  current 
conducting  strips  number  1  and  number  1 1 ,  situated  on  the  edges  of 
the  MEA,  each  of  them  covers  a  smaller  active  area  than  the  other  9 
strips,  collecting  less  currents,  and  the  exact  covering  area  is  hard  to 
calculate  because  it  is  difficult  to  keep  the  edge  of  the  MEA  perfectly 
match  the  edge  of  the  flow  field.  Therefore,  local  currents  number 

I  and  number  11  will  not  be  used  in  the  following  presentation. 
Even  though  the  currents  collected  by  strips  number  1  and  number 

I I  are  not  used  in  the  final  results  analysis,  these  two  strips  were 
still  connected  in  the  measurement  system  to  minimize  extra  cur¬ 
rents  flow  through  the  adjacent  strips.  In  order  to  minimize  edge 
effect  and  focus  on  the  trend  of  the  temperature  distribution,  the 
temperatures  measured  by  T/C  1  and  T/C  10  at  both  ends  will  not 
be  used  in  the  result  presentation,  either. 

Considering  that  this  work  studies  the  current  and  temperature 
distributions  in  steady  state,  average  values  during  a  period  of  5  min 
instead  of  transient  values  at  some  specific  time  will  presented  in 
following  results  for  both  current  distribution  and  temperature  dis¬ 
tribution.  The  5-min-period  starts  only  when  the  local  currents  and 
temperatures  become  stable  and  last  for  a  while. 

Note  that  any  modification  of  PEMFC  for  current  or  temperature 
measurement  may  have  some  influences  on  the  overall  and  local 
performance  of  a  PEMFC.  As  indicated  in  Figs.  1  and  3,  insertion  of 
the  current  distribution  measurement  gasket  will  cause  an  increase 
in  the  gas  channel  depth  which  will  lead  some  change  in  gas  veloc¬ 
ity.  Insertion  of  the  thermocouples  between  catalyst  layer  and  GDL 
in  cathode  could  also  increase  mass  transfer  resistance  in  the  areas 
directly  blocked  by  the  thermocouples,  but  the  thermocouple  is 
so  thin,  and  it  is  inserted  in  under  the  channel,  any  such  effects 
would  be  negligible  since  concentration  losses  in  under  channel 
is  usually  negligible  [29].  Even  if  the  performances  of  the  instru¬ 
mented  cell  would  be  different  from  those  for  an  un-instrumented 
cell,  the  general  relationship  between  local  current  density  and 
local  temperature  would  not  be  affected,  which  is  the  focus  of  this 
study.  Moreover,  compared  with  other  reported  methods,  the  cur¬ 
rent  distribution  measurement  method  used  in  this  work  requires 
minimum  modification  of  a  PEMFC,  and  placing  the  measurement 
gasket  at  the  anode  makes  its  influence  on  the  trend  of  current 
distribution  in  the  in-plane  direction  almost  negligible  according 
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to  a  numerical  study  on  the  influence  of  segmentation  of  current 
collector  by  Natarajan  and  Van  Nguyen  [8]. 

Regarding  the  influence  on  the  temperature  distribution,  which 
is  mainly  due  to  the  lower  thermal  conductivity  of  the  PCB-based 
gasket,  it  can  be  negligible,  too.  Firstly,  it  is  the  temperature  distri¬ 
bution  in  the  in-plane  direction  that  is  studied,  with  structure  and 
material  kept  the  same  for  each  measurement  location.  Secondly, 
it  is  the  steady-state  distribution  that  is  measured.  Therefore,  the 
influence  due  to  lower  thermal  conductivity  of  the  gasket  in  the 
through-plane  direction  would  be  little.  In  addition,  thermocou¬ 
ples  are  placed  between  the  cathode  catalyst  layer  and  the  GDL, 
where  they  are  closest  to  the  electrochemical  reaction  region  gen¬ 
erating  most  heat.  At  this  location,  the  effects  of  other  material  on 
the  measured  temperature  are  further  reduced. 


3.  Results  and  discussion 

With  the  measurement  method  described  above,  current  and 
temperature  distributions  in  the  experimental  PEMFC  under  dif¬ 
ferent  operating  conditions  were  measured  and  analyzed.  In  the 
base  case,  the  flow  rates  for  hydrogen  and  air  are  200  Nml  min-1 
(equivalent  to  stoichiometry  of  1.8  for  average  current  density  at 
1000mA cm-2)  and  1200  Nml  min-1  (equivalent  to  stoichiometry 
of  4.5  for  average  current  density  at  1000mA cm-2),  respectively, 
and  the  anode  and  cathode  humidification  temperatures  as  well 
as  the  cathode  endplate  temperature  are  all  50  °C.  The  flowrates 
chosen  in  base  case  are  quite  large  to  minimize  the  effect  of  insuf¬ 
ficient  reactant.  When  the  focus  is  to  study  the  effects  of  a  specific 
parameter,  only  this  parameter  is  varied  and  all  the  other  parame¬ 
ters  are  kept  unchanged.  For  all  measurements,  both  the  anode  and 
cathode  vents  were  directly  connected  to  atmosphere  to  keep  the 
backpressure  at  1  atm. 


3.1.  Different  air  flow  rates 

Current  and  temperature  distributions  in  the  experimental 
PEMFC  at  different  air  flow  rates  are  plotted  in  Fig.  4.  As  is  shown 
in  Fig.  4(a),  when  the  air  flow  rate  is  very  low  at  60Nmlmin-1, 
equivalent  to  an  average  current  density  of  200  mA cm-2,  the  local 
current  density  decreases  gradually  along  the  direction  from  gas 
inlet  to  gas  outlet.  This  is  caused  by  the  gradual  consumption  of 
oxygen  along  the  flow  direction  in  a  serpentine  flow  field.  The  local 
current  density  becomes  almost  zero  near  the  outlet,  indicating 
oxygen  depletion.  Fig.  4(b)  shows  that  temperature  distribution 
presents  a  similar  trend  of  decreasing  gradually  along  the  flow 
direction  for  air  flow  rate  of  60Nmlmin-1.  As  the  air  flow  rate 
increases,  all  the  local  current  densities  increase  accordingly  due  to 
increased  local  oxygen  concentrations  and  the  current  distribution 
curve  becomes  more  uniform.  This  trend  is  consistent  with  the  pre¬ 
vious  results  on  current  distribution  measurement  in  PEMFC  with 
serpentine  flow  field  [17,49].  It  can  be  seen  from  Fig.  4(b)  that  the 
local  temperatures  also  increase  with  increasing  air  flow  rate  and 
the  temperature  distribution  becomes  more  uniform  similar  to  the 
current  distribution. 

It  is  interesting  to  note  that  when  the  air  flow  rate  increases 
from  900  Nml  min-1  to  1 500  Nml  min-1 ,  as  is  shown  in  Fig.  4(a),  the 
local  current  densities  increase  significantly,  whereas  the  local  tem¬ 
peratures  hardly  change  as  shown  in  Fig.  4(b).  This  phenomenon 
can  be  attributed  to  the  convective  cooling  effect  of  the  air.  Even 
though  heat  generation  rate  is  higher  at  higher  current  density, 
higher  air  flow  rate  also  has  a  higher  heat  transfer  capability.  At 
some  point,  these  two  opposite  effects  on  temperature  at  the  mea¬ 
surement  location  are  equal  and  the  temperature  does  not  change 
significantly. 


Fig.  4.  (a)  Current  distributions  and  (b)  temperature  distributions  at  different  air 
flow  rates. 


3.2.  Different  hydrogen  flow  rates 

Current  and  temperature  distributions  at  different  hydrogen 
flow  rates  were  also  investigated  and  the  results  are  presented  in 
Fig.  5.  At  low  flow  rates  (from  30  Nml  min-1  to  100  Nml  min-1), 
both  local  current  densities  and  temperatures  increase  significantly 
with  hydrogen  flow  rate,  while  at  sufficiently  high  flow  rate  from 
100  Nml  min-1  to  200  Nml  min-1,  both  current  and  temperature 
change  slightly  along  the  flow  direction.  It  is  also  apparent  from 
Fig.  5  that  when  hydrogen  supply  is  insufficient,  it  is  consumed 
quickly  in  the  upstream  and  the  down  stream  becomes  starved  of 
hydrogen.  It  is  very  interesting  to  note  when  the  cell  is  starved  of 
hydrogen,  the  local  current  density  decreases  more  smoothly  than 
the  local  temperature.  This  is  actually  caused  by  the  lateral  cur¬ 
rent  conduction  due  to  the  very  high  current  density  gradient  in 
the  transition  region.  Since  the  local  current  is  collected  after  it 
passes  through  the  GDL  and  the  GDL  has  very  high  lateral  electrical 
conductivity,  when  the  current  generated  under  one  measurement 
strip  is  much  lower  than  its  adjacent  one,  significant  lateral  elec¬ 
trical  conduction  exists  in  the  GDL.  On  the  other  hand,  the  local 
temperature  is  measured  by  thermocouples  inserted  between  the 
catalyst  layer  and  the  GDL,  the  effect  of  lateral  heat  conduction  is 
not  significant.  Therefore,  in  such  situation,  the  local  temperature  is 
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Thermocouple  Numer 

Fig.  5.  (a)  Current  distributions  and  (b)  temperature  distributions  at  different  hydro¬ 
gen  flow  rates. 

actually  a  better  indication  of  local  reaction  rate  than  the  measured 
local  current  density. 

Although  there  are  similarities  between  the  results  at  different 
hydrogen  flow  rates  and  those  at  different  air  flow  rates,  the  differ¬ 
ences  between  them  are  also  apparent.  First,  when  air  flow  rate  is 
low,  the  local  current  density  starts  rather  low  (about  460  mA  cm-2 
for  air  flow  rate  of  60  Nml  min-1,  shown  in  Fig.  4)  and  gradu¬ 
ally  decreases  along  the  flow  direction.  While  when  the  hydrogen 
flow  rate  is  low  at  30  Nml  min-1,  the  local  current  density  starts 
much  higher  (about  800  mA cm-2,  shown  in  Fig.  5)  and  decreases 
much  more  abruptly  along  the  flow  direction.  This  difference  can 
be  attributed  to  the  much  faster  hydrogen  reduction  kinetics  and 
higher  hydrogen  diffusivity. 

Second,  the  local  current  densities  increase  markedly  when 
air  flow  rate  increases  from  900  Nml  min-1  (equivalent  to 
3 100  mA  cm-2)  to  1500  Nml  min-1  as  is  shown  in  Fig.  4(a),  while 
they  increase  very  little  when  hydrogen  flow  rate  increases  from 
100  Nml  min-1  (equivalent  to  900  mA cm-2)  to  200  Nml  min-1 
as  is  shown  in  Fig.  5.  It  suggests  that  hydrogen  flow  rate  of 
100  Nml  min-1  is  already  sufficient  for  the  experimental  PEMFC 
under  such  operating  conditions  due  to  the  faster  transfer  and 
reaction  of  hydrogen  than  oxygen  in  the  air  and  further  increasing 
hydrogen  flow  rate  does  not  improve  the  cell  performance. 


Thermocouple  Numer 

Fig.  6.  (a)  Current  distributions  and  (b)  temperature  distributions  at  different  cell 
operating  voltages  with  sufficient  air  supply. 

Fig.  5(b)  shows  that  the  local  temperatures  with  hydrogen  flow 
rate  of  200  Nml  min-1  are  lower  than  those  with  hydrogen  flow 
rate  of  1 00  Nml  min-1 .  This  is  similar  to  the  case  when  air  flow  rate 
increases  from  900  Nml  min-1  to  1 500  Nml  min-1  discussed  above. 
In  this  case,  increase  in  hydrogen  flow  rate  does  not  leads  to  a  higher 
current  density,  but  the  higher  flow  rate  causes  higher  convective 
heat  transfer,  thus  causes  lower  local  temperature. 

3.3.  Different  cell  voltage 

It  is  well  known  that  operating  voltage  has  significant  effects  on 
both  overall  and  local  current  density  of  PEMFC.  Lower  operating 
voltage  means  higher  overpotential,  which  results  in  higher  current 
density  and  lower  efficiency,  both  causes  higher  heat  generation 
rate.  Therefore,  current  distribution  and  temperature  distribution 
in  the  experimental  PEMFC  under  different  operating  voltages  were 
investigated.  Fig.  6  shows  the  results  with  sufficient  air  supply  of 
1200  Nml  min-1  and  Fig.  7  the  results  with  insufficient  air  supply 
of  only  1 20  Nml  min-1 . 

As  is  shown  in  Fig.  6(a),  all  the  local  current  densities  increase 
significantly  as  the  operating  voltage  decreases  and  the  local  cur¬ 
rent  density  decreases  noticeably  from  inlet  to  outlet  at  low 
operating  voltages,  which  agrees  well  with  previous  results  on  cur¬ 
rent  distribution  measurement  [5,17]. 
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Fig.  7.  (a)  Current  distributions  and  (b)  temperature  distributions  at  different  cell 
operating  voltages  with  insufficient  air  supply. 

Lower  operating  voltage  results  in  higher  local  current  density 
and  higher  local  heat  flux,  and  local  temperatures  should  increase 
accordingly.  This  trend  is  verified  by  the  temperature  distribu¬ 
tion  results  shown  in  Fig.  6(b),  in  which  all  the  local  temperatures 
increase  significantly  with  decreasing  operating  voltages,  similar  to 
the  variation  of  local  current  densities.  Comparison  of  Fig.  6(a)  and 
(b)  shows  that  the  temperature  distribution  trend  generally  corre¬ 
sponds  to  that  of  current  distribution,  indicating  good  correlations 
between  them. 

It  can  be  seen  from  Fig.  7  that  the  results  with  insufficient  air 
supply  show  similar  trend  to  those  in  Fig.  6:  local  current  densi¬ 
ties  increase  with  decreasing  operating  voltage  and  variation  of 
local  temperature  distribution  follows  that  of  current  distribution  h 
very  well.  The  differences  between  Figs.  7  and  6  are  that  local  cur¬ 
rent  densities  and  local  temperatures  in  Fig.  7  are  lower  and  they 
decrease  more  quickly  along  the  flow  direction  due  to  much  lower 
oxygen  concentration  and  lower  gas  flow  velocity. 

It  can  also  be  seen  from  Fig.  6  that  when  the  local  current  den¬ 
sities  are  about  1000  mA  cm-2,  the  local  temperatures  measured 
at  the  cathode  catalyst  layer/GDL  are  about  5°C  higher  than  the 
cathode  endplate  temperature,  TCep ,  which  is  often  taken  as  the 
fuel  cell  operating  temperature  in  many  studies.  When  oxygen  is 
used  as  oxidants  and  PEMFC  operates  at  higher  current  density, 
the  temperature  difference  should  be  even  larger.  This  result  indi- 


Fig.  8.  Correlation  between  current  distribution  and  temperature  distribution  at 
different  air  flow  rates. 

cates  that  caution  must  be  taken  when  using  the  temperature  in 
the  end  plates,  or  the  collector  plate  to  monitor  fuel  cell  operations, 
especially  in  high  current  density  region  and/or  at  high  operating 
temperatures. 

3.4.  Correlation  between  local  current  density  and  temperature 

The  results  presented  in  Figs.  4-7  show  that  the  trends  of  local 
temperature  and  current  density  distributions  are  very  similar.  To 
clearly  show  the  relationship  between  them,  plots  showing  tem¬ 
perature  versus  current  density  are  to  be  produced.  However,  as 
shown  in  Figs.  1  -3,  local  temperatures  are  measured  under  the  flow 
channels  and  local  currents  are  collected  above  flow  field  shoul¬ 
ders,  i.e.  the  locations  for  temperature  and  current  measurements 
are  staggered.  In  order  to  plot  temperature  verses  current  density, 
the  measured  local  temperatures  were  interpolated  linearly  to  the 
positions  of  local  current  densities.  Additionally,  instead  of  local 
temperature,  the  temperature  rise,  AT,  defined  as  the  difference 
between  the  measured  local  temperature  and  TCep  (50  °C  in  this 
work),  is  used  in  the  presentations. 

Figs.  8  and  9  show  the  relationships  between  the  local  temper¬ 
ature  rise  and  local  current  density  at  the  same  cell  voltage  and  at 


Fig.  9.  Correlation  between  current  distribution  and  temperature  distribution  at 
different  hydrogen  flow  rates. 
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Fig.  10.  Correlation  between  current  distribution  and  temperature  distribution  at 
different  cell  voltages  with  sufficient  air  supply. 

different  air  and  hydrogen  flow  rates,  respectively.  It  can  be  clearly 
seen  that  the  local  temperature  rise  increases  almost  linearly  with 
current  density  for  every  specific  gas  flow  rate  when  the  cell  volt¬ 
age  is  kept  constant.  It  can  be  seen  also  that  the  slope  of  curves  at 
different  flow  rates  are  almost  the  same. 

The  linear  correlation  can  be  explained  with  basic  thermody¬ 
namic  analysis  on  heat  generation  in  PEMFC.  It  is  known  that  local 
electric  energy  output  (Pei)  and  local  heat  generation  (q)  in  PEMFC 
can  be  simply  expressed  as 

Pei=jU  (1) 

Q  —  ji^thermo  ~  U)  (2) 

where  j  represents  the  local  current  density,  U  the  cell  operating 
voltage,  and  Ethermo  the  thermo  neutral  potential,  which  is  about 
1.48  V  for  hydrogen  PEMFC  at  standard  temperature  and  pressure 
(STP)  based  on  the  high  heating  value  of  hydrogen. 

Eqs.  (1)  and  (2)  suggest  that  if  the  cell  operating  voltage  is 
kept  constant  everywhere,  which  is  true  for  the  case  shown  in 
Figs.  8  and  9  with  a  fixed  operating  voltage  of  0.3  V,  the  local  heat 
generation  is  proportional  to  the  local  current  density.  It  shows 
that  the  correlation  results  can  be  explained  and  agree  well  with 
the  basic  thermodynamic  analysis. 

It  can  be  noticed  by  careful  observation  of  Fig.  8  that  slopes  of 
the  correlation  curve  decrease  as  air  flow  rate  increases.  It  can  be 
attributed  to  the  increased  convective  cooling  effects  of  higher  air 
flow  rate. 

Figs.  10  and  11  show  the  correlation  results  at  different  cell 
operating  voltages  with  sufficient  air  supply  ( 1 200  Nml  min-1 )  and 
insufficient  air  supply  ( 1 20  Nml  min-1 ),  respectively.  It  can  be  seen 
that  slopes  of  the  correlation  curve  increases  as  the  operating  volt¬ 
age  decreases.  The  trend  is  more  apparent  in  Fig.  1 1  than  that  in 
Fig.  10  due  to  the  low  air  flow  rate  of  120  Nml  min-1,  where  the 
convective  cooling  effects  are  not  as  significant.  This  trend  can  also 
be  explained  by  the  basic  thermodynamic  analysis  shown  in  Eq.  (2). 
According  to  Eq.  (2),  at  the  same  local  current  density,  the  lower  the 
operating  voltage  is,  the  higher  is  the  heat  generation  rate,  thus  the 
slope  of  the  temperature  versus  current  density  curves  increases 
with  decreasing  cell  voltage. 

Such  good  correlations  between  local  temperature  and  current 
density  could  provide  foundations  for  novel  in  situ  diagnostic  tech¬ 
niques  in  fuel  cells  since  local  temperature  is  usually  much  easier 
to  measure  and  it  could  be  used  to  determine  local  current  density 
[32], 


Fig.  11.  Correlation  between  current  distribution  and  temperature  distribution  at 
different  cell  voltages  with  insufficient  air  supply. 


4.  Conclusions 

Combined  with  the  current  distribution  measurement  gasket 
technique  and  high-precision  thin  thermocouples  placed  at  the 
cathode  catalyst  layer/GDL  interface,  a  technique  for  simultane¬ 
ous  measurements  of  local  current  density  and  temperature  in 
an  PEMFC  is  developed.  With  this  technique,  simultaneous  mea¬ 
surements  are  realized  in  a  commercially  available  experimental 
PEMFC.  From  the  results  of  simultaneous  measurements  under  dif¬ 
ferent  operating  conditions,  the  following  conclusions  can  be  made : 

•  There  exist  good  correlations  between  local  temperature  and 
local  current  density. 

•  The  correlations  between  local  temperature  and  local  current 
density  agree  very  well  with  basic  thermodynamic  analysis. 

•  The  slope  of  the  curves  for  local  temperature  rise  versus  local 
current  density  increases  with  decreasing  cell  operating  voltage. 

•  The  effect  of  convective  cooling  on  local  temperature  in  PEMFC  is 
significant  at  high  reactant  flow  rates. 

•  Significant  temperature  difference  exists  between  cathode  cata¬ 
lyst  layer  and  the  endplate  plate. 

•  Simultaneous  measurement  of  current  and  temperature  distri¬ 
butions  can  provide  more  detailed  information  than  separate 
measurements. 
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